Hole Dispersion and Symmetry of the Superconducting Order Parameter for 
Underdoped CuC>2 Bilayers and 3D Antiferromagnets. 
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correlations. 

The purpose of this paper is (i) the study of the proper- 
ties of holes and (ii) the extension of previous analysis for 
the symmetry of the superconducting order parameter to 
the case of lightly doped bilayers and 3D antiferromag- 
nets. The effective Hamiltonian introduced in previous 
publications for one plane H is defined on the 2D square 
lattice as 
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We calculate the dispersion of a hole in C11O2 bilayers and 3D antiferromagnets, using the self- 
consistent Born approximation (SCBA) to the t — J model. Superconductivity and the symmetry 
of its order parameter are studied introducing a nearest-neighbor density-density attraction induced 
by short range antiferromagnetic (AF) fluctuations, as described in recent studies for the single 
layer cuprates. The well-known pairing in the d x i_ y i channel observed for one plane remains robust 
when three dimensional interactions are turned on. In bilayers, as the exchange along the direction 
perpendicular to the planes grows, eventually a transition to a "s-wave" state is observed with 
an order parameter mixture of d iz i_ r i and s x 2 +y 2 +z 2 . For an isotropic 3D antiferromagnet the 
d x 2_ y 2 and d 3z 2_ r 2 channels are degenerate. Our results are compared with other predictions in the 
literature. 

Pacs Numbers: 74.20.-z, 74.20.Mn, 74.25.Dw 
I. INTRODUCTION 



Since the discovery of high temperature superconduc- 
tors, the origin of the pairing mechanism has been contro- 
versial. Numerous studies suggest that the normal state 
properties are anomalous and, as a consequence, above 
the critical temperature, T c , the strongly correlated char- 
acter of the conduction electrons in the cuprates should 
be taken into account in any realistic study. In particu- 
lar, the presence of short-range antiferromagnetic fluctu- 
ations have been invoked by several groups as the origin 
of the deviations from a conventional Fermi liquid behav- 
ior observed in transport measurements above T c . Theo- 
ries based on antiferromagnetic fluctuations quite clearly 
predict a superconducting condensate in the d x 2_ y 2 chan- 
nel. EH These studies are supported by experimental 
results for the bilayer cuprate YBa2Cu307_^, which are 
compatible with a highly anisotropic pairing state, likely 
a d x 2_ y 2 singlet. || Photoemission results for Bi-2212, 
which also have two CuOi layers per unit cell, support 
this scenario, |I| and thus superconductivity induced by 
AF correlations is a strong candidate to explain the pair- 
ing mechanism of the cuprates. 

The authors and collaborators have recently explored 
these ideas using simple models for quasiparticles (q.p.) 
considered as holes strongly dressed by AF fluctuations, 
and interacting through a nearest-neighbor (NN) density- 
density attraction, as suggested by the study of the two- 
holes bound-state of the t — J model. @] The presence of 
a large accumulation of weight in the density of states, 
induced by the small bandwidth of the q.p. dispersion, 
enhances T c beyond what would be naively expected from 
the strength of the hole-hole attraction. These ideas were 
originally referred to as the "antiferromagnetic van Hove" 
(AFVH) scenario, but it is important to remark that the 
boost in T c does not arise exclusively from a van Hove 
singularity but mainly from the "flat" regions observed 
in the q.p. dispersion, effect which is induced by AF 
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where Ck Q denotes a destruction operator of a q.p. with 
dispersion extracted from accurate numerical and analyt- 
ical studies of one hole in the undoped t — J model given 
by eAF(k)/eV = 0.165 cos k x cos k y + 0.0435(cos2fc x + 
cos2fcj,), where t — OAeV and J/t = 0.3 are assumed. 
The quasiparticles move within the same sublattice 
to avoid distorting the AF background (a = A,B de- 
notes the sublattice). This is correct in the limit of a 
small number of holes in a long-range ordered AF sys- 
tem, and it is expected to be a good approximation even 
with AF short-range order, as long as the AF correla- 
tions are strong. The dispersion e^(k) reflects a re- 
markable feature of the cuprates, observed in photoemis- 
sion experiments, namely the presence of flat bands near 
(0, 7r) (it, 0) in the standard square lattice notation. |||7]] 
The parameter V is the intersublattice density-density 
attractive coupling strength between holes also discussed 
in previous literature, || with ri\ the number operator. 
The rest of the notation is standard. 

The interaction in Eq.(l) should be considered as the 
dominant piece of a more general and extended AF- 
induced effective potential between holes which in the 
presence of long-range order corresponds to a sharp 8- 
function centered at Q = (jr, n) in momentum space, 
and it acquires a width as the AF correlation length 
decreases. Even with £af as small as a couple of lat- 
tice spacings, it has been shown that the NN interaction 
used in Eq.(l) remains robust. S Equivalently, the NN 
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hole-hole attraction can be considered as arising from the 
"minimization of broken AF links" argument ^| in the 
large J/t limit. 

It is expected that Hamiltonian Eq.(l), which is cer- 
tainly a very simplified version of the low energy behav- 
ior of the t — J model, nevertheless captures the im- 
portant features of holes moving in AF backgrounds at 
low doping, and temperatures smaller than J where the 
quasiparticles are dominant. This is a regime difficult 
to study with numerical methods directly applied to the 
t — J model. The analysis of Eq.(l) using standard BCS 
techniques, Exact Diagonalization approaches, and the 
Eliashberg equations, has been fairly successful in repro- 
ducing some experimental features of the cuprates. [^]|]] 
For example, in this model superconductivity appears in 
the d x 2_ y 2 channel, as it seems to occur in experiments, 
with a critical temperature T c ~ 1QQK . In addition, the 
concept of "optimal doping" appears naturally due to the 
robust peak in the density of states (DOS) produced by 
the dispersion 6af0<-), which is crossed by the chemical 
potential as the hole density grows. In this respect, the 
theory has features very similar to those discussed be- 
fore in "van Hove" theories for the cuprates. [|l0| Indeed, 
the quasiparticle lifetime is linear with energy at optimal 
doping. However, the flat regions in the dispersion and 
the associated accumulation of weight in the DOS are 
produced by strong correlations and thus they are sta- 
ble against small perturbations, like impurities and extra 
couplings, effects that usually destroy weak logarithmic 
van Hove singularities. The specific heat jump, 2A/T C 
and other important BCS ratios are in good agreement 
with the experimental data. Q 

Although certainly more work is needed to show that 
these "real-space" pairing ideas presented in previous lit- 
erature are a strong candidate to describe the cuprates, 
its quantitative success led us here to study geometries 
and couplings beyond those of the single layer with some 
confidence. Our goal is to report the trends observed 
when systems different from a single layer cuprate are 
analyzed. Here, the hole spectrum in the bilayer and 3D 
antifcrromagnets is calculated with the SCBA which was 
shown to reproduce accurately numerical results for the 
2D t — J model, [[ll) and it was successfully applied to the 
2D t — t 1 — J model to compare its predictions with the 
photoemission spectra of Sr2Cu02Cl2- |l2| Using the for- 
malism previously described, a hole-hole attraction will 
be introduced producing a superconducting state. Here, 
the symmetry of this superconducting state will be an- 
alyzed and compare with other predictions for the same 
bilayer and 3D systems. 



II. SELF CONSISTENT BORN 
APPROXIMATION. 

The Hamiltonian for spins and holes used in this sec- 
tion is defined as, 

h = - X] '"E^i"^'' + h - c ^ 

ex. ia - 

+ ^Ja^[(S?Sf +ea -inin i+e J 

a i 

+^ a (S?Sr +ea +SrS+ +e J], (2) 

where i denote sites of a bilayer or simple 3D cubic clus- 
ter, e a is the unit vector in the a-direction (a = x, y, z), 
t a and J a correspond to the NN hopping amplitude and 
exchange coupling, respectively, in the direction a, t a p is 
the next-nearest-neighbors (NNN) hopping in the direc- 
tions defined by e a ± ep, the parameter ^ Q ,(e[0, 1]) rep- 
resents a possible exchange anisotropy added for com- 
pleteness, and the rest of the notation is standard. In 
the terms with NNN interactions, the summation is done 
with the condition a / ft i.e. we only consider NNN 
hopping along the diagonals of the plaquettes. A gen- 
eralization to include NNN hopping at distance of two 
lattice spacings along the main axes is straightforward. 
The constraint of no double occupancy is implemented 
in the kinetic energy by means of the standard definition 
for hole operators c\„ — Ci<j(l — rii_ CT ), with n- la = c ia Ci a . 

Let us extend the SCBA analysis which has been ap- 
plied before to single planes (TIJ to the bilayer and 3D 
Hamiltonian Eq.(2). Following steps that by now are 
standard, first we redefine operators according to Sj — > 

S"? — > — *^T' c j ct — * c j-<7! where j denotes sites on the 
B sublattice. Spins on the A sublattice remain intact. 
With this procedure an AF state becomes ferromagnetic. 
Then, the usual linearized Holstein-Primakoff transfor- 
mation is used. This transformation is defined in terms 
of Bose operators a\ and a[, according to w ai^/2S, 
Sr w a\V2S, and S? = S - a\cn. To handle the 
hole sector, the spin-charge decomposition q-f = h\ and 

C U = ^H^t i s introduced, where hi, h\ are operators 
corresponding to spinless holes. Finally, using periodic 
boundary conditions and after long but straightforward 
algebra, we can rewrite the resulting Hamiltonian in mo- 
mentum space as 

k q 
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where the Bose operators fr(b') and q(a^) are connected 
by a Bogoliubov transformation, JTl] TV is the num- 
ber of sites in the lattice, ek = 4(t xy cos k x cos k y + 
txz cos kj^ cosfc 2 + t yz cos ky cos k z ) is the bare spinlcss 
hole dispersion (which cancels in the absence of an ex- 
plicit NNN hopping in the Hamiltonian), and Mk iq = 
M qAc-q + VqPk is the hole-magnon vertex. The magnon 
dispersion is given by w q = wo^q, where ujq = zSJ, 
J = ^ J2 a J°<> & = \ ls the spin of the electrons in the 
original t — J model language, d is the number of dimen- 
sions (which will be taken equal to 3 in the rest of the 
paper), z = 2d is the coordination number, and we have 
used units where H = 1. The other quantities are defined 
as follows: 

/3k = 2 \_. ta cos k a , 



'1 + V a 
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The ground state energy Eq has the form: 

E = -dJN[S 2 +S(l-iE "q) + ll- 



(4) 



(5) 



Care must be taken in order to apply Eqs.(3) and (4) to 
the bilayer. The small size in the z-direction (one lattice 
spacing) does not allow us to apply periodic boundary 
conditions (BC) to simulate bulk effects, as we do in the 
x- and y-directions. The z-axis BC should be open. How- 
ever, for a bilayer open BC are equivalent to periodic BC 
but with parameters J z , t z , t xz , and t yz reduced in half 
from their actual values. Remember also that k z can 
only be or tt. Thus, for a bilayer assumed symmetric 
with respect to the x <-> y interchange in the planes, and 
spin symmetric ($ Q = 1), two branches of the magnon 
spectrum are found, namely 
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where q is a 2D vector, J = J x = J y is the in-plane 
coupling, J± = J z is the interlayer coupling, and r q = 
2(cos<7 x + cosg y ). fig}] . 

With this information, we can find the quasiparticlc 
hole spectrum in the AF background which is defined by 
the position of the peak with the lowest binding energy 
(assuming its intensity Z is finite) in the spectral function 
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e k — E(k, uj + iS) 
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where the self-energy £(k, u>) is calculated self- 
consistently from 
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which corresponds to the sum of magnon rainbow di- 



agrams, as described in previous literature. |11| Eq.(8) 
allow us to consider, at least in part, the dressing of the 
hole by spin excitations, and it has been shown that the 
predictions arising from this approach are in excellent 
agreement with Exact Diagonalization numerical results 
for 2D finite clusters. jll| Thus, we are confident that 
the rainbow approximation may maintain its accuracy 
when applied to bilayer and 3D systems where numerical 
results are difficult to obtain. The solution of Eq.(8) is 
obtained numerically on finite but large clusters. 

While the SCBA formalism described here was setup 
for arbitrary values of the many couplings appearing in 
Hamiltonian Eq.(2), in the results reported in the fol- 
lowing sections, we have neglected NNN hoppings in all 
directions and set the magnetic anisotropy parameter $ a 
to l(spin isotropic system). We also assumed xy sym- 
metry in the planes, i.e. orthorhombic distortions are 
not considered here. Thus, in the rest of the paper we 
will use the following notations: J and t will refer to 
the in-plane exchange and hopping parameters, while Jj_ 
and t± correspond to the inter-plane ones. The gener- 
ality of Eqs.(2-8) would allow the interested reader to 
obtain SCBA results for arbitrary situations where spin 
or lattice anisotropies are important, or including NNN 
hopping amplitudes. 

Regarding the numerical values of the couplings, we 
have introduced relations between exchange and hopping 
amplitudes to simplify the multiparameter analysis. In 
particular, we have here assumed = (^") 2 , which 
arises from the t — J model when derived from the one 
band Hubbard model in strong coupling. In addition, we 
have fixed t — 0.4 eV throughout the paper. Thus, as 
parameters we only have J/t and J±/J. 

The actual value of the inter-plane exchange is contro- 
versial. The two-magnon peak in Raman spectra exper- 
iments for YBa2Cu30g+x is consistent with an in-plane 
J ~ 125 meV, with a negligible interplane exchange. H] 
Recent measurements of the nuclear quadrupole and nu- 
clear magnetic resonances on Cu-sites in YBa2Cu30i5 
(alternating 1-2-3 and 1-2-4 structures) suggest a ratio 
^f- ^0.3. However, a different value ~0.55 was 
found from infrared transmission and reflection measure- 
ments, after applying linear spin wave theory. fl6| ] A hop- 
ping amplitude ratio ^f- = 0.4 was reported on the basis 
of first-principles linear density approximation calcula- 
tions, pH corresponding to ^j- ~0.16 with our conven- 
tion. In view of these discrepancies, in most of the results 
below the perpendicular exchange is consider a free pa- 
rameter. 
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III. RESULTS 



A. Bilayer. 

Let us first discuss the one hole spectrum in the AF 
bilayer system using the exchange ratio — (^j-) 2 = 
0.16 suggested by band structure calculations, t = 0.4eV 
and J = 0.125 eV, as found in some experiments (J/t = 
0.3125). The self-consistent calculation was performed on 
a 16 x 16 x 2 cluster. In Fig.l the spectrum is presented. 
It has two branches that correspond to the two possible 
momenta in the z-direction i.e. and n. These branches 
obey the relation e + (k) = e_(k + Q), where Q = (it, n) 
and k is the 2D momentum (e±(k) here denote the even 
and odd branches of the spectrum). This relation is in- 
duced by the presence of antiferromagnetic long-range 
order in the system, which effectively reduces the size of 
the Brillouin zone. The minima of the branches are close 
to, but not exactly at, k = (tt/2,tt/2) i.e. where the 
minimum of the one-layer dispersion resides. This small 
effect is observed from fits of the data using trigonomet- 
ric functions and it is barely noticeable in Fig.l, but as 
we increase Jj_ the splitting of the minima away from 
k = (7r/2,7r/2) becomes larger, and at J± — J it is 
clearly visible. The bilayer one-hole spectrum can be 
accurately fit by NN, NNN and next to NNN hoppings 
amplitudes, namely e±(k) = ±0.01(cosfc2; +cosfc y ) + 
0.099 cos k x cos k y + 0.033(cos2fc x + cos2fc a ) (eV) (see 
Fig.l). Note that each individual branch contains a NN 
amplitude which is nonzero. This is necessary since in- 
dividually the even and odd branches are not invariant 
under a momentum shift in Q. It should also be remarked 
that near (ir, 0) there are "flat" dispersion regions, as ob- 
served in previous studies of the t — J model. |6[0] It has 
been argued that AF correlations induce these features 
in models of correlated electrons even when the ^af is 
only a couple of lattice spacings. jl8| Fig.l shows that 
these features also survive the introduction of a realistic 
bilayer coupling. 

Our SCBA calculations indicate that the antiferromag- 
netic bilayer has quasiparticle excitations since a large 
peak appears at the bottom of the hole spectral function 
A(k, ui). Typical spectral functions for the high symme- 
try points are presented in Fig. 2. The quasiparticle peak 
carries a substantial percentage of the spectral weight 
in agreement with the calculations for a single layer. ||] 
Note also that the closer the momentum to the bottom 
of the band, the stronger the quasiparticle peak. We ob- 
served a similar behavior even for interlayer exchanges as 
large as ij- — 3. 

The evolution of the Fermi Surface (FS) as the SCBA 
hole dispersion is populated in the rigid band picture 
is shown in Fig. 3. At very low hole density it starts 
with hole pockets around (7r/2,7r/2) (Fig. 3(a)) for the 
even and odd branches. The pockets are longer along 
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(0,0) (71/2,71/2) (71,0) (7T,0) (0,0) 

FIG. 1. SCBA spectrum of one hole in an AF bilayer for 
parameters if = (if) 2 = 0.16, t = 0.4eV and J = 0.125 
eV, on a 16 x 16 x 2 cluster. Open circles belong to the bonding 
branch e+, while stars correspond to the antibonding branch 
e_. Results between (ir/2,n/2) and (7r,0) (not shown) are 
trivially related by symmetry to others presented in the same 
figure. The solid (dot-dashed) line is a fit using trigonometric 
functions for the bonding (antibonding) branch. 
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co(eV) 

FIG. 2. SCBA spectral functions of a hole in an AF bi- 
layer at high symmetry points for the same values of param- 
eters used in Fig.l, and with an artificial broadening S = O.li 
for the 5-functions. The momenta (k x ,k y ,k z ) are indicated. 
k z — corresponds to the bonding branch, while k z = 7r to 
the antibonding. 
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the (0, tt) — (7T, 0) direction than along the main diag- 
onal in the Brillouin zone as observed in similar stud- 
ies for the single layer problem. jl9|,[3| As these pockets 
grow in size, the Fermi level eventually hits the saddle- 
points and the FS changes its topology (Fig. 3(b) and 
(c)), becoming a large FS when the chemical potential 
is above the saddle points (Fig. 3(d)). It should be re- 
marked that after the hole pockets disappear, the FS 
acquire quasi-nesting features, which may lead to an en- 
hancement of certain susceptibilities. The energy scale 
of the rapid change from hole pockets to a large FS is 
about 400K. This implies that a strong temperature de- 
pendence should be expected in transport measurements 
of a doped AF bilayer, similar to those observed in a sin- 
gle AF layer. Jl^] Actually the overall FS evolution ob- 
tained from the rigid band filling of the bilayer dispersion 
for ^j- = {^y) 2 — 0.16 is qualitatively similar to results 
found for just one CuC"2 plane using the same approxima- 
tions, Jl^H with the main difference being the presence 
of two branches. Identifying similar features in experi- 
ments could give us information about the strength of 
the bilayer coupling and antiferromagnetic correlations. 





FIG. 3. Evolution with doping of the Fermi Surface ob- 
tained by a rigid band filling of the one hole SCBA spectrum 
presented in Fig.l. The density grows moving from (a) to (d). 
The thick (thin) line corresponds to the bonding (antibond- 
ing) branch. 



B. 3D Cubic Lattice 

In this subsection, the properties of a hole in a 3D 
antiferromagnetic environment will be analyzed. The 
SCBA quasiparticle hole dispersion for the case of an 



isotropic cubic lattice (i.e. 



j 



it) = 1) i s shown in 



Fig. 4. The numerical self-consistent calculation was car- 
ried out on an 8 x 8 x 8 cluster with parameters chosen 



which is natural based on the results reported for the 
one and two coupled planes. Our best fit of the nu- 
merical data corresponds to the following dispersion: 
e(k) = 0.082(cos k x cos k y +cos k y cos fc 2 +cos k x cos k z ) + 
0.022(cos2/c ;E + cos2k y + cos 2^) (eV). As in the case of 
the one layer system, holes tend to move within the same 
sublattice to avoid distorting the AF background. 
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FIG. 4. SCBA spectrum of one hole in an isotropic 3D AF 
with J/t = 0.4. The momentum in the z-direction is in- 
dicated. The lines are fits to the numerical data: the solid 
line corresponds to k z = 0, dashed to k z — tt/2, and dot- 
ted-dashed to k z = tt. Points between (tt/2, tt/2) and (n,0) 
(not shown) are related by symmetry to others presented in 
the figure. 

It is interesting to analyze the dependence of the quasi- 
particle dispersion total bandwidth, W, with the ratio j 
for the isotropic 3D AF. Here W is defined as the differ- 
ence in energy between the state with the highest energy, 
typically at (0,0,0) and (tt,tt,tt), and the state with the 
lowest energy at (tt/2, tt/2, tt/2). Fig. 5 shows W as the 
ratio j varies from 0.1 to 1. Results for a single layer 
on an 8 x 8 cluster are also shown for comparison. The 
overall behavior of the bandwidths is qualitatively the 
same in both cases, with W for a 3D AF being slightly 
larger. In both cases W/t grows approximately linearly 
for small J/t, and reaches saturation around j ~ 0.6. 
In the linear regime, we find W ~ 2.5 J showing that the 
characteristic energy scale of the dispersion is J rather 
than t, a well-known result obtained before in the con- 
text of holes in 2D antiferromagnets. It is interesting 
to observe that the increase in the dimensionality of the 
problem does not change dramatically the bandwidth of 
holes. 

As for planes and bilayers, the SCBA calculations show 
a strong quasiparticle peak in the hole spectral function 
A(k, uj) of an isotropic 3D AF. In Fig. 6 results are shown 
at particular high symmetry points. The parameters are 
the same as those used to calculate the dispersion of 
Fig. 4. We checked that the relative weight of the quasi- 



as t = 0.4eV and J = 0.16eV, which gives j = 0.4. 
The minimum of the dispersion is at (tt/2, tt/2, tt/2), particle peak increases with the ratio j. The quasiparti- 
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FIG. 5. Open circles represent the bandwidth W vs J/t 
obtained from the SCBA quasiparticle hole spectrum for an 
isotropic 3D AF. Filled diamonds are the bandwidth of an 
isotropic 2D AF (evaluated on an 8 x 8 cluster) shown for 
comparison. Lines are drawn to guide the eye. 

cle weight Z is similar in 2D and 3D systems, at least at 
the bottom of the hole spectrum. 
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FIG. 6. Spectral functions of a hole in an isotropic 3D AF 
at high symmetry points for the same values of parameters 
used in Fig.4 with a broadening 5 = O.lt for the 5-functions. 
The 3D momenta k is indicated. 



IV. SUPERCONDUCTIVITY. 

If there is a source of attraction between the quasipar- 
ticles described in the previous sections, then the ground 
state of the system could become superconducting. In 
the underdoped regime, AF correlations are strong and, 
thus, it is natural to study pairing mediated by AF fluc- 
tuations assuming that the dispersion at half-filling is not 
drastically distorted when a finite but small hole density 



is studied. Since £af is likely of the order of only a cou- 
ple of lattice spacing at realistic dopings, a real — space 
approach to pairing is more suitable than a picture where 
extended states (magnons) are interchanged between car- 
riers (k — space approach). While these two variations of 
AF-mediated pairing are quantitatively different, never- 
theless it has been shown that they lead to the same sym- 
metry for the superconducting order parameter (SCOP). 
H Thus, independent of the discussion of real-space vs 
k-space approaches, we believe that it is safe to use the 
real-space pairing ideas of Ref. || to study the symmetry 
of the SCOP for the case of a bilayer and 3D AF systems. 

As explained before, the simplified picture introduced 
in previous literature || is to use the spectrum of a hole, 
calculated numerically or using the SCBA, as the disper- 
sion of the quasiparticles which interact through a NN 
attraction regulated by a parameter V (of the order of 
J). This interaction is motivated by AF-induccd pairing 
since it can be shown that a potential V(q) in momen- 
tum space maximized at q = Q and smeared by a finite 
£af, induces a real-space potential which is dominated 
by a NN attraction. Previous estimations suggest a value 
V = 0.6 J from the study of the two-holes bound state on 
a single layer. Q Here, the same value will be employed 
for the in-plane NN interaction, while for the interaction 
in the z-direction V± = 0.6 J± will be used. In the absence 
of accurate Exact Diagonalization or Monte Carlo results 
for the dispersion of holes away from half-filling, a rigid- 
band filling of the SCBA quasiparticle spectra obtained 
at half-filling will be assumed in the BCS gap equation. 
Previous studies for single layer systems have shown that 
this approach provides a SCOP with a symmetry consis- 
tent with other more traditional methods, and in addition 
the actual values of T c and its density dependence are in 
qualitative agreement with experiments. Q Thus, it is 
natural to employ the same approach to analyze the SC 
properties of doped bilayers and 3D antiferromagnets. 



A. Superconductivity in a bilayer. 

Using the bilayer SCBA hole dispersion and the hole- 
hole NN attraction described before, we have studied su- 
perconductivity within the BCS formalism. To calcu- 
late the SCBA hole dispersion, we fixed j — 0.3125 
and allow for the interlayer coupling to take the values 
^f- =0.16, 0.3, 0.5, 1, 2 and 3. Solving numerically the 
gap equation on large clusters, we observed that for the 
first three ratios ^j- the symmetry of the superconduct- 
ing order parameter (SCOP) is d x 2_ y 2 for all the hole 
densities where superconductivity exists. This is the re- 
gion analytically connected to the physics of one layer, 
where a similar condensate was found in previous papers. 
In Fig. 7, we show T c vs the quasiparticle hole density for 
a bilayer with ^j- = 1 . This dependence is also shown for 
the case of a single layer with the same J/ 1, for compar- 
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ison. Both curves exhibit the "optimal doping" behavior 
and high values of T c of the order of 100K. The differ- 
ence in T c between a bilayer and a single layer is small 
for practical purposes, except for the region of high den- 
sities where a bilayer changes the symmetry of the super- 
conducting condensate into an "s-wave" state to be dis- 
cussed below (note that a quasiparticle weight Z smaller 
than one implies the existence of incoherent weight in 
the spectral function at energies larger than those of the 
q.p. band. Then, "high density" in this context does 
not imply a very small electronic density (n)). Similar 



conclusions were obtained in Ref. |20| where an Exact 
Diagonalization study of coupled layers was carried out 
to calculate the superconducting pairing correlations. No 
important enhancement of pairing correlations changing 
from a single t — J layer to a bilayer was reported. pC|| 
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FIG. 7. Superconducting critical temperature T c (solid 
line) vs quasiparticle hole density for an AF bilayer with 
J/t = 0.3125 and J±/J = 1 (see text for details). The 
dashed line is T c vs quasiparticle hole density for a 2D AF 
plane with J/t = 0.3125. Both lines are calculated within 
the framework of the real-space pairing approach and the BCS 
gap equation. For bilayers, at high q.p. densities (x ~ 55%) 
there is a transition to an "s-wave" state. 

According to the classification of Liu, Levin and Maly 
pl| the d x 2_ y 2 state that we observed is labeled (d,d), 
which means that a gap with the symmetry d x 2_ y 2 opens 
in both the bonding and antibonding Fermi surfaces (i.e. 
the Fermi surfaces corresponding to k z — and ir shown 
in Fig.l). This d-wave channel dominates up to values 
close to Jj_ ~ J. At the ratio ^j- =1 the symmetry 
of the SC order parameter changes only at high densi- 
ties (i.e. around 55% filling of the quasiparticle band as 
shown in Fig. 7) to a mixture of d 3z 2_ r 2 and s x 2 +y 2 +z 2, 
with real amplitudes of the same sign. s x 2 +y 2 +z 2 is the 
3D analog of the extended s-wave function. The rela- 
tive weight of the g? 3z 2 _ r 2 component is dominant in the 
bilayer case. In general, the tetragonal point symmetry 
group consists of 10 irreducible representations among 



which only five correspond to singlet pairing. The func- 
tions 3z 2 — r 2 , x 2 + y 2 + z 2 and also a constant function 
have the same transformation properties with respect to 
the tetragonal symmetry operations and thus belong to 
the same representation of the tetragonal lattice sym- 
metry group. Then, they mix together and they form the 
"s-wave" representation with respect to rotations in the 
planes, while d x 2_ y 2 belongs to and d xy to for 
a tetragonal system, p^] Note that in an orthorhombic 
lattice G?3 Z 2_ r 2, s x 2 +y 2 +z 2 as well as d x 2_ y 2 correspond to 
the same representation T\. 

A less trivial issue is how the symmetry evolves from 
a mixture of d^ z 2_ r 2 and s x 2 +y 2 +z 2 to the so-called d z 
state [[23) (the d z state is a singlet along the link in the 
z-direction of the bilayer, i.e. it is "s-wave" with respect 
to 7r/2 rotations of the planes). This transition occurs 
as the ratio ^j- grows at a fixed density. The relative 
weight of s x 2 +y 2 +z 2 in the mixture increases until both 
functions become equally weighted with the same signs 
of the weight coefficients in the limit of J± ^> J. In this 
situation the mixture is exactly equivalent to the d z sym- 
metry. In our calculations a coupling ^j- = 2 is enough 
to have almost d z symmetry for all densities, with a sub- 
stantial increase of T c up to 150K. The mixture of d3 Z 2_ r 2 
and s x 2 +y 2 +z 2 can also be viewed as the mixture of d z 
and the extended in-plane s-wave s x 2 +y 2, or in terms of 
odd-even gaps A±(k) cx ±a + fl(cosk x +cosk y ), where 
k is the in-plane momentum of the quasiparticle (the even 
and odd gaps differ by a constant corresponding to the 
inter-plane pairing amplitude). In a similar fashion, as 
the ratio ^j- increases, the relative weight of the in-plane 
extended s-wave decreases. 
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FIG. 8. Schematic representation of the strong V coupling 
limit of the two quasiparticle problem for the case of a bilayer. 
The dashed and dotted lines represent the effective hopping 
amplitudes contained in the SCBA dispersion that we label 
t\ and t-2, respectively. For more details see text. 

There is a simple intuitive way to understand the evo- 
lution with J±/J of the symmetry of the superconduct- 
ing order parameter found here. j^] Consider the limit of 
Eq.(l) applied to a bilayer where the NN in-plane and 
inter-plane attractions dominate over the kinetic energy. 
In this limit, let us analyze the two quasiparticle problem 
and its bound state. The quasiparticles will be roughly 
at a distance of one lattice spacing (tight bound state). 
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In the reference frame of one particle, the other can oc- 
cupy any of the five NN sites (four in the same plane and 
one in the other plane). These sites are actually linked 
by the dispersion in the SCBA since the five sites be- 
long to the same sublattice. In other words, the wave 
function of one quasiparticle orbiting around the other 
can be found in this limit from a simple 5-site 1-particle 
problem as shown schematically in Fig. 8, with the ra- 
tio of hopping amplitudes t\ and £2 simulating the effect 
of Jj_/J- Solving this trivial problem we find that for 
J±/J < 1 (or £2/^1 < 1) the lowest energy is obtained 
for a state with amplitudes c, — c, c, — c, (c > 0) at sites 
1,2,3,4,5, respectively, which is the d x 2_ y 2 state. On 
the other hand, for J±/J > 1 (or t 2 /ti > 1) the ground 
state amplitudes have the same value (with its sign) in- 
plane, and a different value and sign out of the plane, 
as expected from the discussion above. Note also that 
these simple calculations illustrate the fact that a NN 
attraction and on-site repulsion is not enough to have a 
tendency to form d-wave pairs in a dilute gas. In the 
present case this potential must be supplemented by a 
particular hole dispersion that allows for the movement 
of carriers within the same sublattice. 

B. Comparison with other theories 

In the literature, several types of interactions between 
quasiparticles in a bilayer have been proposed. Of partic- 
ular interest for our discussion are those based on phe- 
nomenological magnetic susceptibilities in the presence 
of dynamical short-range AF order, both in and between 
the planes. Recently it has been suggested |24|] that the 
SCOP in a bilayer can be of a particular s-wave type 
named (s, — s) with opposite signs in the odd and even 
branches of the spectrum of carriers. The assumption 
made in Ref. |Q is that the in-plane magnetic suscepti- 
bility X|| (q) has the opposite sign, but the same absolute 
value, as the inter-plane susceptibility XJ-(q)- in other 
words, X||(q) = -X±(q) = Xo(q), where xo(q) is the 
phenomenological susceptibility of a doped 2D antifer- 
romagnet peaked at q = Q = (tt, tt) introduced to fit 
NMR data. |25) This assumption needs further justifica- 
tion since the analysis of experimental results provides an 
estimation 5 meV < J± < 20 meV E(| for the interplane 
coupling, while J w 120 meV. Thus, it is somewhat risky 
to assume that the spin correlations between planes have 
the same strength as in the planes. Then, it seems more 
natural to employ the RPA approximation for the suscep- 
tibility m given by xi*(q) = Xo^i^TJx, where x±(q) 
are the even and odd branches of the bilayer susceptibil- 
ity, respectively. In this case, X||(q) an d XJ-(q) are quite 
different in absolute value, i.e. |x||(q)l ^ IX-i(q)l m 
order to satisfy the applicability of the RPA approxima- 
tion. We have checked by a direct calculation, that this 
approach provides d x 2_ y 2 symmetry for the SCOP on 



both branches of the spectrum i.e. the so-called (d, d) 
state. This is in agreement with Ref. pl| where a simi- 
lar phenomenological in-plane susceptibility derived from 
the neutron scattering data was used to construct the ef- 
fective potential with an arbitrary relation between the 
in-plane and inter-plane coefficients. The above men- 
tioned RPA model is certainly relevant to the more phe- 
nomenological approach used in this paper if one consid- 
ers the expansion of xo(q) in real space, assuming that 
the effective coupling between quasiparticles is propor- 
tional to xo- For correlation lengths £ ~ 2 — 3 lattice 
spacings, the expansion consists of a strong on-site re- 
pulsion U, and also a relatively strong NN attraction V 
in-plane, with a ratio j^y ~ 3. The NNN interaction 
is repulsive and about one order of magnitude weaker 
than U. Other terms are negligibly small. Then, the 
phenomenological interaction used here and in previous 
literature || is very closely related to the RPA interac- 
tion in the realistic regime Jj_ <C J. To show this result, 
first let us expand the susceptibility as 

= -if TT r ~ * o(q) T J ^o(q)- (9) 
Xo (<l)±J± 

Using a 3D notation one can write xi^lz) — Xo(q) — 
cos(q z ) JiXo(q)' f° r the case of a bilayer with q z — or tt. 
The effective in-plane potential must contain an strong 
on-site repulsion. Consequently, Xo(q) can be written as 
Xo(q) = u + /(q), where u is a constant proportional to 
the on-site repulsion U, and the function /(q) satisfies 
the condition ^ q /(q) = 0. Hence, the susceptibility 
can be rewritten as %(q, q z , ) ~ u + /(q) — cos(g 2 ) J±[u + 
/(q)] 2 . Assuming that the ratio is small (i.e. the 
on-site repulsion is strong), we finally get: 

X(q, + /(q) - J±u cos(q z ) - 2J±u cos(g z )/(q). 

(10) 

The first three terms are accounted for in the calcula- 
tion discussed in this paper both for a bilayer and a 3D 
antiferromagnet since, as it was argued above, the largest 
contribution to /(q) comes from the NN attraction. The 
last term in Eq.(10) is relatively small compared to the 
other three, but it is assumed to have a significant con- 
tribution in the models of Refs. explaining the 
difference between their results (i.e. (s, — s)) and ours 
(i.e. mixture of d 3z 2_ r 2 and s x -2 +y 2 +z i). Nevertheless, 
both cases are "s-wave" in-plane and thus there is no 
qualitative difference in their transformation properties 
(they belong to the same representation). In addition, 
since J±/J <C 1, both are unlikely to be stable in the 
cuprates. 

C. Superconductivity in the isotropic 3D AF. 

To calculate T c and the symmetry of the SCOP in the 
case of a three dimensional isotropic AF, we here use the 
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SCBA hole dispersion obtained at J/t = 0.4 and again 
the BCS gap equation. As expected, |27],|2£| when ^f- = 
1, the two channels d x 2_ y 2 and d 3z 2_ r 2 are degenerate, 
i.e. T c is the same for both and varies with the hole 
density as it is shown in Fig. 9. The critical temperature 
in 3D is slightly smaller than for planes and bilayers. The 
"optimal doping" behavior is also present in 3D, but it 
is less pronounced than in 2D. 
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FIG. 9. Superconducting critical temperature T c vs quasi- 
particle hole density for an isotropic 3D AF with J/t — 0.4, 
using the BCS gap equation and Hamiltonian Eq(l). 

When ^j- ^ 1, the system acquires tetragonal sym- 
metry and the channels d x 2_ y 2 and d 3z 2_ r 2 are no longer 
degenerate. In fact, if ^4- < 1, the gap possesses d x 2_ y 2 
character, while if ^j- > 1, it becomes again a mixture of 
d 3z 2_ r 2 and s x 2 +y 2 +z 2 , as it was mentioned in the case of 
a bilayer. Again, as the ratio increases, the symmetry 
of the SCOP becomes d z . All these results can also be 
obtained following a discussion similar to that of Fig. 8 
but now using six sites instead of five. 



V. CONCLUSIONS 

In this paper, we have investigated the properties of 
bilayer and 3D antiferromagnets specially regarding the 
behavior of holes in such environments and the even- 
tual formation of a superconducting state upon doping. 
The main assumption is that AF correlations remain 
of at least a few lattice spacings in range which, as a 
first approximation, allow us to use the dispersion of 
holes at half-filling and a density-density NN attraction 
between quasiparticles also generated by antiferromag- 
netism. Within the BCS gap equation superconductiv- 
ity is investigated. A small coupling in the z-direction 
does not alter qualitatively the d x 2_ y 2 superconducting 
state found for a pure single layer plane. Only for large 
values of the exchange in the direction perpendicular to 



the planes is that a competing "s-wave" superconducting 
state becomes dominant. 



VI. ACKNOWLEDGMENTS 

A. N. acknowledges the support by the National High 
Magnetic Field Laboratory, Tallahassee, Florida. E. D. 
is supported by the NSF grant DMR-9520776. We thank 
MARTECH for additional support. 



[1] 
[2] 
[3] 

[4] 
[5] 
[6] 
[7] 

[8] 
[9] 
[10] 

[11] 

[12] 

[13] 
[14] 

[15] 
[16] 
[17] 

[18] 
[19] 
[20] 



D. J. Scalapino, Phys. Rep. 250, 331 (1995), and refer- 
ences therein. 

E. Dagotto, Rev. Mod. Phys. 66, 763 (1994); and refer- 
ences therein. 

D. A. Wollman, D. J. Van Harlingen, W. C. Lee, D. M. 
Ginsberg, and A. J. Leggett, Phys. Rev. Lett. 71, 2134 

(1993) ; C. C. Tsuei, et al., Phys. Rev. Lett. 73, 593 

(1994) . 

Z.-X. Shen et al, Science, 267, 343 (1995); J. Ma et al., 
Phys. Rev. B 51, 3832 (1995); and references therein. 

E. Dagotto, A. Nazarenko and A. Moreo, Phys. Rev. 
Lett. 74, 310 (1995). 

E. Dagotto, A. Nazarenko and M. Boninsegni, Phys. Rev. 
Lett. 73, 728 (1994). 

N. Bulut, D. J. Scalapino and S. R. White, Phys. Rev. B 
50, 7215 (1994); R. Preuss, W. Hanke, and W. von der 
Linden, Phys. Rev. Lett. 75, 1344 (1995); M. Langer, J. 
Schmalian, S. Grabowski, and K. H. Bennemann, Phys. 
Rev. Lett. 75, 4508 (1995). 

A. Nazarenko, S. Haas, J. Riera, A. Moreo, and E. 
Dagotto, preprint. 

A. Nazarenko, A. Moreo, J. Riera, and E. Dagotto , to 
appear in Phys. Rev. B (RC). 

C. C. Tsuei et al., Phys. Rev. Lett. 65, 2724 (1990); R. 
S. Markiewicz, J. Phys. Condens. Matt. 2, 6223 (1990); 
A. A. Abrikosov et al., Physica C 214, 73 (1993). 
S. Schmitt-Rink, CM. Varma and A.E. Ruckenstein, 
Phys. Rev. Lett. 60, 2793 (1988); F. Marsiglio, et al., 
Phys. Rev. B 43, 10882 (1991); G. Martinez and P. 
Horsch, Phys. Rev. B 44, 317 (1991); Z. Liu and E. 
Manousakis, Phys. Rev. B 45, 2425 (1992). 
A. Nazarenko, K. J. E. Vos, S. Haas and E. Dagotto, R. 
J. Gooding, Phys. Rev. B 51, 8676 (1995). 
N.E. Bonesteel, Phys. Rev. B 47, 11302 (1993). 
K. B. Lyons et al, Phys. Rev. Lett. 60, 732 (1988); S. 
Sugai et al, Phys. Rev. B 42, 1045 (1990); 
R. Stern et al, Phys. Rev. B 50, 426 (1994). 



M. Griininger et al, preprint cond-mat/9501065 



O. K. Andersen, O. Jepsen, A. I. Liechtenstein, and I. I. 

Mazin, Phys. Rev. B 49, 4145 (1994). 

A. Moreo et al., Phys. Rev. B 51, 12045 (1995). 

S. Trugman, Phys. Rev. Lett 65, 500 (1990). 

J. A. Riera, Phys. Rev. B 49, 3629 (1994) 



9 



[21] D. Z. Liu, K. Levin, and J. Maly, Phys. Rev. B 51, 
(1995). 

[22] M. Sigrist and T. M. Rice, Z. Phys. B. 68, 9 (1987). 
[23] E. Dagotto, "Computational Approaches in Condensed 

Matter Physics", Eds. S. Miyashita, M. Imada, and H. 

Takayama, Springer- Verlag Proceed, in Physics Vol.70; 

N. Bulut, D. J. Scalapino and R. T. Scalettar, Phys. 

Rev. B 45, 5577 (1992); E. Dagotto, J. Riera, and D. 

Scalapino, Phys. Rev. B 45 5744 (1992); R. Eder, Y. 

Ohta, and S. Maekawa, Phys. Rev. B 52, 7708 (1995) 
[24] A. I. Liechtenstein, I. I. Mazin, and O. K. Andersen, 

Phys. Rev. Lett. 74, 2303 (1995). See also R. J. Radtke, 

A. I. Liechtenstein, V. M. Yakovenko, and S. Das Sarma, 

preprint; S. Grabowski, J. Schmalian, M. Langer, and K. 

H. Bennemann, preprint. 
[25] A. J. Millis, H. Monien, and D. Pines Phys. Rev. B 42, 

167 (1990). 

[26] A. J. Millis and H. Monien, preprint [rxmd-mat/950608g . 
[27] K. Miyake, S. Schmitt-Rink, and C. M. Varma, Phys. 

Rev. B 34, 6554 (1986) 
[28] D. J. Scalapino, E. Loh, Jr., and J. E. Hirsch, Phys. Rev. 

B 34, 8190 (1986) 



10 



